Monte Carlo Study of the Spin-1 Baxter- Wu Model 
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The two-dimensional spin-1 Baxter- Wu model is studied by using Monte Carlo simulations. The 
standard single-spin-flip Metropolis algorithm is used to generate the configurations from which 
the order parameter, specific heat and magnetic susceptibility are measured. The finite-size scaling 
procedure is employed in order to get the critical behavior. The extensive simulations shown that 
the critical exponents are different from those of the spin-1/2 model suggesting that the spin-1 model 
is in a different universality class. 

PACS numbers: 64.60Kw, 64.60Cn, 64.60Fr 



INTRODUCTION 

The Baxter- Wu model is a system of spins defined on a 
two-dimensional triangular lattice with the classical spin 
variables Si taking only integer values. It was first intro- 
duced by Wood and Griffiths I'] as a model which does 
not exhibit invariance by a global inversion of all spins. 
The system is described by the Hamiltonian 



n 



<ijk> 



(1) 



where the coupling constant J is positive and the sum 
is over all triangles made up of nearest-neighbor sites on 
the triangular lattice. For the spin-1/2 model, where 
Si = ±1, the exact solution obained by Baxter and 
Wu gives fcsTe/J = 2/ln(l V2) and a = = §0. 
The system has also been studied with quenched impu- 
rities by Monte Carlo yl and Monte Carlo renormaliza- 
tion group approaches Conformal invariance stud- 
ies [sl [Q] have shown that the pure spin-1/2 Baxter- Wu 
and the four-state Potts models have the same opera- 
tor content and are in the same universality class. More 
recently, the short time critical dynamics has been in- 
vestigated through the relaxation of the order parameter 
at the critical temperature by Monte Carlo simulations 
0. On the other hand, for spin values greater or equal 
to one there are neither exact solutions nor even much 
approximate approaches. It is the purpose of this work 
to study the model above for the spin-1 case by using 
Monte Carlo simulations, where the variables Si take the 
values Si = —1, 0, 1. 

Monte Carlo methods 0, 01 form the largest and most 
important class of numerical methods used for solv- 
ing statistical physics problems. The basic idea behind 
Monte Carlo simulation is to simulate the random ther- 
mal fluctuation of the system from state to state over the 
course of an experiment. Performing a high-precision 
finite-size scaling analysis using standard Monte Carlo 
techniques is very difficult due to constraints on the 
available computer resources. The introduction of his- 



togram techniques to extract the maximum information 
from Monte Carlo simulation data at a single temper- 
ature enhances the potential resolution of Monte Carlo 
methods substantially 10. 11]. In this sense, we apply 
the histogram techniques together with the Metropolis 
simulation algorithm in order to investigate the thermal 
behavior of the spin-1 Baxter- Wu model defined by Eq. 

by considering the specific heat, order parameter and 
magnetic susceptibility. Our main interest is to obtain, 
through a finite-size scaling analysis, the phase transi- 
tion temperature as well as the critical exponents of the 
model. 

In the next section we present the thermodynamic 
quantities and the details of the simulations. In section 
III we discuss the results and in section IV we summarize 
our conclusions. 



SIMULATION BACKGROUND 

The simulations have been carried out by using the 
single-spin-flip Metropolis algorithm. In the course of the 
simulations we considered triangular lattices with linear 
dimensions L x L and fully periodic boundary conditions 
for system sizes of length 18 < L < 108. Due to the fact 
that the system has, in addition to the ferromagnetic 
phase (with all spins up), three different ferrimagnetic 
phases with three different sublattices (one sublattice up 
and spins on the other two sublattices down) the allowed 
values of L are always a multiple of 3. In this way, all 
ground states of the infinite lattice would fit on any finite 
lattice. Following equilibration (which comprised 6 x 10"* 
MCS) runs comprising up to 5 x iC^ MCS (Monte Carlo 
steps per spin) were performed. Histogram reweighting 
|lOl im and finite-size scaling techniques were used to 
precisely locate the second-order phase transition. Re- 
garding the histograms, great care has been taken in or- 
der to assure the reliabily of the extrapolated results for 
all lattice sizes. 

The thermodynamic quantities we measured in our 



2 



simulations are the order parameter, defined as the root 
mean square average of the magnetization of the three 
sublattices 



niB 



mc 



(2) 



where mA, tub and mc are the magnetizations per spin 
of the different sublattices, the order parameter suscep- 
tibility defined as 



(3) 



where (3 — l/ksT (with ks the Boltzmann constant and 
(...) means an average over the generated Monte Carlo 
configurations), and the specific heat 



(4) 



where {E) is the mean value of the energy. 

According to finite-size scaling theory the critical tem- 
perature scales as 



Tr, 



(5) 



where A is a constant, Tc is the critical temperature of 
the infinite system, and is the effective transition tem- 
perature for the lattice of linear size L. This effective 
temperature can be given by the position of the maxi- 
mum of any of the following quantities: the temperature 
derivative of m. In to or Inm^, the order parameter sus- 
ceptibility or the specific heat. The above temperatures 
are given in units of J /ks- An independent estimate of 
V, however, can be made through the evaluation of the 
maximum logarithmic derivative of any power of the or- 
der parameter to" since one has 



dU_ 



where a is a constant and U is either In to or Ini 



(6) 



(or, 



in general. In to"). In addition, the specific heat and the 
magnetic susceptibility scale, at the transition tempera- 
ture, as 



(7) 



where a and 7 are the critical exponents of the specific 
heat and susceptibility, respectively. From Eqs. (|2I7|I one 
can obtain the critical temperature and critical exponents 
of the model. 



RESULTS 

The independent evaluation of the critical exponent 
V, as obtained from Eq. ^ without any consideration 
of the critical temperature Tc, is shown in Fig. ^ for 



the maximum derivative of the logarithm of m and m? 
(although other powers of to can also be used). From 
both cases one has v = 0.617(4), which is different from 
V = 0.666 for the spin-1/2 model. 

With V determined quite accurately we proceed to es- 
timate the position of Tc- As discussed in the previous 
section, the location of the maxima of the various thermo- 
dynamic derivatives, namely the maximum of the specific 
heat, susceptibility, and the derivatives of to and In to and 
InTO^, provide estimates for the transition temperature 
which scale with system size as Eq. (jSJ. A plot of these 
estimates is given in Fig. |21 The results from the lin- 
ear fit are listed in Table HI One can note that they are 
indeed quite close to each other and a final estimate is 
Tc = 1.6607(3). 

The logarithm of the maximum value of the specific 
heat and order parameter susceptibility as a function of 
the logarithm of L is shown in Fig. O From these data 
one has a = 0.692(8) and 7 = 1.13(1). 




FIG. 1: Logarithm of the maximum values of the derivatives 
of In m and In m? as a function of the logarithm of the size L. 
The straight lines, with slopes corresponding to v — 0.617(3) 
in both cases, show the asymptotic behavior of the fits. The 
errors are smaller than the symbol sizes. 



TABLE I: Estimated critical temperatures from different ther- 
modynamic quantities according to the linear fit shown in Fig. 



Quantity 



c 

X 

( dm \ 
\ dT 'max 
' d In m A 

dT J^max 
d In m'' 



1.6607(1) 
1.6605(1) 
1.6606(1) 
1.6609(1) 
1.6610(2) 
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1.690 - 



1.680 




0.000 



0.002 



0.004^.,,^, 0.006 



0.008 



0.010 



FIG. 2: Size dependence of the effective critical tempera- 
tures (in units of J/ks) estimated from several thermody- 
namic quantities. The lines are fits to Eq. @ with v = 0.6i7 
obtained from Fig. Qand the intercepts are given in Table in 
The errors are smaller than the symbol sizes. 
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FIG. 3: Logarithm of the maximum values of the specific 
heat C and order parameter susceptibility x as a function of 
the logarithm of L. The straight lines are fits to Eqs. ||7J with 
a/v = 1.121(6) and 7/1/ = 1.829(9). The errors are smaller 
than the symbol sizes. 



CONCLUSIONS 

It is clear, from the quality of the above resuhs, that 
a well defined second order phase transition takes place 
in the model at Tc = 1.6607(3) with critical exponents 
y = 0.617(3), a = 0.692(6) and 7 = 1.13(1), which are 
indeed different from the spin-1/2 case, namely v = 2/3, 
a ~ 2/3 and 7 = 7/6. This means that this three spin in- 
teraction model has exponents which depend on the spin 
value. It is worth saying that the present model can also 
have an extra interaction with a crystal field of the form 



D^-Si^. This is a generalization in the direction of the 
so-called Blume-Capel model 0|. What we have done 
here is studied the special case D = Q. However, in the 
limit D —00 one recovers the spin-1/2 model. From 
the present results we then expect that along the second- 
order transition line for different values of D one has a 
line with varying critical exponents. In addition, as we 
have shown, a second-order phase transition takes place 
at = in contrast with the conjecture that the spin-1 
Baxter- Wu model is critical only in the limit Z) — > — 00 
[1^. Some preliminary results, agreeing with the picture 
of a line of second-order phase transition with varying ex- 
ponents and the presence of a multicritical point, for the 
present system with crystal field interaction, have already 
been obtained from conformal invariance with finite-size 
scaling theory and the mean field renormalization group 
approach [lj|. Work in this direction using Monte Carlo 
simulations is now in progress. 
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